Two-color QCD in a strong magnetic field: The role of the Polyakov loop 
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We study two-color QCD in an external magnetic backround at finite temperature using the 
Polyakov-loop extended two-flavor two-color NJL model. At T = 0, the chiral condensate is cal- 
culated and it is found to increase as a function of the magnetic field B. In the chiral limit the 
deconfinement transition lies below the chiral transition for nonzero magnetic fields B. At the phys- 
ical point, the two transitions seem to coincide for field strengths up to \qB\ « 5m, whereafter they 
split. The splitting between the two increases as a function of B in both the chiral limit and at the 
physical point. In the range from zero magnetic field and \qB\ = 20m^, the transition temperature 
for the chiral transition increases by approximately 35 MeV, while the transition temperature for 
deconfinement is essentially constant. 

PACS numbers: 



I. INTRODUCTION 



The behavior of hadronic matter at finite tempera- 
ture and density in strong external magnetic fields has 
received a lot of attention for many years, see for ex- 
ample Ref. [l[ for a very recent review. The problem of 
strongly interacting matter in a strong magnetic back- 
ground arises in various contexts. For example, magne- 
tars, which are a certain type of neutron stars, have very 
strong magnetic fields of the order of 10 10 Tesla (T) [2J. 
Some of the properties of stars such as the mass-radius 
relation are determined by the equation of state. The 
determination of the bulk properties of a Fermi gas in an 
external magnetic field is therefore important for the un- 
derstanding of these compact stellar objects. Similarly, 
large magnetic fields, up to the order of eB ~ iq 14-16 
T, where e is the electric charge of the pion. are being 
generated in noncentral heavy-collisions at the Relativis- 
tic Heavy-Ion Collider (RHIC) and the Large-Hadron 
Collider (LHC) [111. The presence of strong magnetic 
fields may be observed in these experiments via the chi- 
ral magnetic effect. This effect is basically the separation 
of charge in a magnetic background due to the existence 
of topologically nontrivial configurations in the decon- 
fined phase of QCD p|. Finally, we mention that strong 
magnetic fields of the order 10 14 — 10 19 T may have been 
present in the early universe during the strong and elec- 
troweak phase transitions [U,0] ■ The presence of a strong 
magnetic field at the electroweak phase transition may 
have implications for baryogenesis, i. e. for the genera- 
tion of the baryon asymmetry in the universe [all . 

Chiral symmetry of the QCD Lagrangian and the 
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spontaneous breaking of this symmetry in the QCD vac- 
uum is an essential feature of the strong interactions. 
At T — 0, it is expected that a constant magnetic 
background enhances chiral symmetry breaking if it is 
present already at B — or that it induces chiral sym- 
metry breaking if the symmetry is intact at B = 0. This 
phenomenon is called magnetic catalysis and has been 
discussed in Refs. 

[TM3 

in the context of the Nambu- 
Jona-Lasinio (NJL) model, chiral perturbation theory, 
and QED (note however the recent paper [l8| where the 
authors argue that effects from the neutral mesons might 
show magnetic inhibition if B is strong enough). The 
basic mechanism is that neutral quark-antiquarks pair 
minimize their energy by both aligning their magnetic 
moments along the direction of the magnetic field [l|. 
Magnetic catalysis was recently demonstrated on the lat- 
tice by Braguta et al 11911 i n three-color quenched QCD 
as well as by Bali et al [20J in the context of three-color 
and 1+1+1 flavor QCD. The results of [20|, which are 
for physical quark masses and extrapolated to the con- 
tinuum limit, are reproduced very well up to magnetic 
fields of the order eB ~ 0.1 (GeV) 2 in chiral perturba- 
tion theory [H, [l6[ and up to eB ~ 0.25 (GeV) 2 using 
the Polyakov-loop extended NJL model (PNJL) [2lj . 

Magnetic catalysis at T = gives rise to the expecta- 
tion that the critical temperature T c for the chiral transi- 
tion is an increasing function of the magnetic field B. In- 
deed, c/) 4 -theory (22[, chiral perturbation the ory [23l |24| 
(However, see also Ref . I25l0 . the NJL model |261.[27|. the 
PNJL modeljp. Ill [29^ and the quark-meson (QM) 
model [13, l30l432| all predict this behavior. Further- 
more, the PNJL model also predicts a modest split of 
approximately 2% between the chiral transition and the 
deconfinement transition, except for Ref. (33|. In this 
case the split is of the order 10% and is due to the 
effects of dimension 8 operators. However, bag-model 
calculation [34J , the Polyakov-loop extended QM calcu- 



lation 35J, and the large- N c calculation [36] all predict 
a decreasing critical temperature as a function of B. In 
Refs. [UHll, it is probably related to their treatment of 
vacuum fluctuations and related renormalization issues. 

Turning to lattice simulations, t he p icture seems to 
be complicated as well. In Refs. [33, |38[, the lattice sim- 
ulations indicate that the chiral critical temperature is 
increasing as a function of the magnetic field. In this 
case, the bare quark masses used correspond to a pion 
mass in the range = 200 — 480 MeV, i. e. a very 
heavy pion. These results have been confirmed by Bali 
et al [39]. However, for light quark masses that cor- 
respond to the physical pion mass of = 140 MeV, 
the critical temperature is a decreasing function of the 
magnetic field B [40| . The basic mechanism seems to be 
that the magnetic catalysis at T = turns into inverse 
magnetic catalysis fill |42| for temperatures around the 
critical temperature T c . The results suggest that the 
critical temperature is a nontrivial function of the quark 
masses. 

Two-color QCD is interesting for a number of rea- 
sons. For example, in contrast to three-color QCD, one 
can perform lattice simulations at finite baryon chem- 
ical potential fis- The reason is that due to the spe- 
cial properties of the gauge group SU(2), the infamous 
sign problem is absent and thus importance sampling 
techniques can be used. Moreover, the physics at finite 
baryon chemical potential is very different from its three- 
color counterpart: Again due to the properties of the 
gauge group, two quarks can form a color singlet and so 
diquarks are found in the spectrum of the chirally broken 
phase. The diquarks are bosons and finite baryon chem- 
ical potential is then the physics of relativistic bosons 
and their condensation at low temperature. In the chi- 
ral limit, the Lagrangian of two-color two-flavor QCD 
has an 5C/(4) symmetry. Since this group is isomorphic 
to 5*0(6), chiral symmetry breaking can be cast into the 
form 5*0(6) — > 50(5). The Goldstone modes are there- 
fore contained in a single five-plet with the usual three 
pions, a diquark and an antidiquark as well. Various 
aspects of the phase diagram of two-color QCD can be 
found e. g. in Refs. |43l - [52l ]. 

The problem of two-color QCD in a strong mag- 
netic background was first investigated on the lattice by 
Buidovidovich, Chernodub, Luschevskaya, and Polikar- 
pov [53L |54| in the quenched approximation. Magnetic 
catalysis at T = has been verified and in in the chi- 
ral limit, the chiral condensate grows linearly for small 
values of B. This behavior is in qualitative agreement 
with chiral perturbation theory. Later, lattice simula- 
tions have been carried out with dynamical fermions by 
Ilgenfritz, Kalinowski, Miiller-Preussker, Petersson, and 
Schreiber [55j. Their results seem to indicate that the 
condensate grows linearly with B in the chiral limit at 



T = 0. They also found that for all temperatures and 
fixed bare quark mass, the chiral condensate grows with 
the magnetic field. This implies that the critical tem- 
perature is an increasing function of the magnetic field. 

In the present paper, we use the PNJL model to study 
two-color QCD in a constant magnetic background B 
at finite temperature and zero baryon chemical poten- 
tial. The article is organized as follows. In Sec. II, we 
briefly discuss the PNJL model in a magnetic field and 
the thermodynamic potential. In Sec. Ill, we present 
our numerical results and in Sec. IV, we summarize and 
conclude. 



II. PNJL MODEL AND THERMODYNAMIC 
POTENTIAL 

In this section, we briefly discuss the two-flavor two- 
color PNJL model. The Euclidean Lagrangian can be 
written as 



£ — Cq + C\ + £ 2 , 
where the various terms are 



(1) 



£ = #7 m £>m - m ]V> ' ( 2 ) 
£1 - d [W) 2 + {^15^) 2 + {^f + ferVO 2 

+ [^a 2 T 2 ^\ 2 + |^ 7 5<T 2 T 2 Vf] , (3) 
£2 = G 2 [{^) 2 - (ipij5'4') 2 - (TpTip) 2 + (ipijuTi/j) 2 

-\W^T 2 ^\ 2 + \Wl5<J2T2^\ 2 , (4) 

where the quark field ip is an isospin doublet 

*-(;)• < 5 > 

The covariant derivative is = — iqA^ — ioiA % ^ 
where is the gauge field associated with £7(1) elec- 
tromagnetism and A 1 ^ is associated with SU{2) color. 
The covariant derivative is diagonal in flavor space, 
q = diag(2/3, -l/3)|e|. a t {i = 1,2,3) are the Pauli 
matrices acting in color space, while Ti are the Pauli 
matrices acting in flavor space, mo is the mass ma- 
trix which is diagonal in flavor space and contains the 
bare quark masses m u and m^. In the following we take 
m u = md- Moreover, tp denotes the charge conjugate 

of the Dirac spinor, ip c = Cip , where C = iy 2 ^ . G\ 
and G 2 are coupling constants. The interacting part £1 
is invariant under global 11(A) — 5C/(4) x ^(l)^ trans- 
formations while the £2 is invariant under global 5t/(4). 
One sometimes writes G\ = (1 — a)G and G 2 = aG 
and so the parameter a determines the degree of ^7(1)^ 
breaking. In the following we choose a—\. 



2 



We next introduce the collective or auxilliary fields 

a = -2Grpip ,7Tj = -2G$i'j 5 T i '>p , pi = -2Gipnip , 

A = -2G^7 5 T 2 cr 2 i/. c ,A 5 = -2G^iT 2 a 2 ip c , (6) 

where er, iii, pi, A, and A5 have the quantum num- 
bers of a scalar isoscalar, pseudoscalar isovector, scalar 
isovector, scalar diquark, and psuedoscalar diquarks, re- 
spectively. The Lagrangian (JlJ can then be written com- 
pactly as 

£ = [i^D^ - ra - er - i^n^ - r iy Oj] ip 

A*ip c i"/ 5 T 2 a 2 ip + H. c. + A* 5 ip c 'ir 2 a 2 ijj + H. c. 



now bilinear in the quark fields and we integrate them 
out exactly by performing a Gaussian integral. This 
gives rise to an effective action for the composite fields. 
In the mean-field approximation, we neglect the fluctua- 
tions of the composite fields and the fermionic functional 
determinant reduces to 



2 
1 

AG 



IAJ 2 



(?) 



If we use the equation of motion for er, 7r$ t;, A, and A5 
to eliminate the auxilliary fields from the Lagrangian ([7]), 
we obtain the original Lagrangian (fT]). 

In pure gauge theory, the Polyakov loop $, which is 
the trace of the Wilson line l{x) = e i ^o dr ai A\(x,r) ^ j e 
$ = -^-Tre^o 3 dT<TiA 4 (x,T) ^ - 1S an orc j er parameter for de- 
confinement [56| . Under the center symmetry Zn c , it 
transforms as $ -> e 27rn / N ^ j where n = 0, 1, ...N c — 1. 
For N c = 2, this is simply a change of sign and in two- 
color QCD the Polyakov is purely real. At low temper- 
ature, i. e. in the confined phase, we have $ « and in 
the deconfined phase, we have 

In the PNJL model, a constant background temporal 
gauge field Aa is introduced via the covariant derivative 
in Eq. (P) [13, HI]. In Polyakov gauge, the background 
field is diagonal in color space and, A\ = OS 1 ' 3 , where 9 
is real. The Wilson line can then be written as l(x) = 
diag(e ll3e ,e~ tl3e ) and the order parameter $ reduces to 



$ = cos(/30) 



(8) 



In order to allow for a chiral condensate, we introduce a 
nonzero expectation value for the field er 1 



-2G(ipip) + a , 



(9) 



where a is a quantum fluctuating fields with vanishing 
expectation values. To simplify the notation, we intro- 
duce the quantity M which is defined by 



M = m -2G{ipip) 



(10) 



Note that the expectation value M is assumed spacetime 
independent in the remainder of this paper. Eq. (JT]) is 



1 Since we consider the case of zero quark chemical potential, the 
other collective fields have zero expectation value. 



a 



quark 



-AN n 



{M-m f 
AG 
d 3 p 



(2ir) 3 ' Ep 



Tlog 



1 



(11) 



where E p = ^Jp 2 + M 2 and E± = E p ± Q. Note that 
the integral involving E p is ultraviolet divergent and re- 
quires regularization. We will return to this issue below. 
Summing over ±, we can write 



quark 



(M-TO ) 2 

AG 



AN r 



d 3 p 
(2tt)3 



E p + Tlog [l + 2<S>e- 0E ? +e 



-2@E P 



(12) 



The interpretation of Eq. (Tl"2"]) is now as follows. For 
$ 0, we have confinement and thus a thermal 
part proportional to Tlogfl + e~ 2l3Ep ], which corre- 
sponds to an excitation of energy 2E p , i. e. a bound 
state. Similarly, for $ w 1, the thermal part is 
Tlog [1 + 2$e- /3B " + e" 2 ^] = 2Tlog[l+e- /3B "] which 
is the thermal contribution from two degrees of freedom 
each with energy E p , i. e. two deconfined quarks. 

The complete thermodynamic potential 51 is given by 
the sum of the contributions from the quarks, f2 qU ark 
in Eq. (111[) and a contribution from the gluons, f2 ga ugo: 
where |48| 



augc 



= -bT [24$ 2 e - /3a + log(l - $ 2 )] , (13) 



where a and b are constants. This form is motivated 
by the lattice strong-coupling expansion [62j . In the 
pure gauge theory, we can find an explicit expression for 

the value of the Polyakov loop, |$| = yjl — ^e' 30 and 

so T c = a log 24. $ goes to zero in a continuous man- 
ner showing that the phase transition is second order in 
agreement with universality arguments [561 ]. 

We next consider this system in a constant magnetic 
field B along the z-axis. We do this by using the 
covariant derivative D,, 



= 5a2XiB and Ai = 5 1 < 3 S 



<9 M - iqa^ - i<JiA 1 ^ where 
^40. Note that the SU(4) 
symmetry of the Lagrangian is broken in an external 
magnetic field due to the different electric charges of the 
u and d quarks. The remaining symmetry is a U{1)a 
symmetry which corresponds to a rotation of the u and 
d quarks by opposite angles [15j. The chiral condensate 
breaks this Abelian symmetry and it gives rise to a single 
Goldstone boson, namely the neutral pion. 
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The energy eigenvalues of the Dirac equations are in 
this case given by 



^p 2 z +M 2 + {2m+l~s)\q f B\ , (14) 



where M is the mass of the quark, s is the spin of the 
quark with electric charge qj and m denotes the mth 
Landau level. In Eq. (|12p . dispersion relation E p is now 
changed to E m and the three-dimensional integral be- 
comes a one-dimensional integral and a sum of Landau 
levels m. For a quark with charge qf, we then make the 
replacements 



p 2 -> vl 



d 3 p 

(27T) 5 



WfB\ 
2tt 



(2m 

E 



l + s)\q f B\ 

dp z 
2ir ' 



(15) 
(16) 



where the sum m is over Landau levels and where the 
prefactor takes into account the degeneracy of the 
Landau levels. The divergent term in Eq. (fT2j) is denoted 
by Slq„^ rk and now becomes 



div 
quark 



= -N r 



dp z 
2tt 



kfB\ 



f,m,s 



M 2 



{2m-l + s)\q f B\ . (17) 



The integral over p z as well as the sum over Landau 



levels m in Eq. (fTT)) are divergent. We will use zeta- 
function regularization and dimensional regularization 
to regulate the divergences. The integral is now gener- 
alized to d — 1 — 2e dimensions using the formula 



d d p 
{2ir) d 



M% = - 



; 7E A 2 

47T 



_1 2 > hj.I 

(47T) 



(18) 



where Mg 



M 2 + (2m + 1 - s)\q f B\ and A is the 
renormalization scale in the MS renormalization scheme. 
This yields 



O div 

J 'quark 



-Nr. 



e7EA2V F(-l + £ ) V«M— 



-in 



Air 2 

f,m,s 



(19) 



For each flavor /, the sum over m and s can be written 

as 



2-2r 
B 



{2q f Bf 



C(-l + £,£/) 



1 



(20) 



where C(a,x) = Y^ 00 , -, — r^rr is the Hurwitz zeta func- 

tion and Xf = 2 |^b | • Expanding the Hurwitz zeta func- 
tion in powers of e, we obtain 



div 



16tt 2 ^ 



A 2 



2\q f B\ 



2(q f B) 



M l 



1 

- + 1 

e 



8(g / i?) 2 C (1 ' 0) (-l,x / )-2| g/ i?|Af 2 log: 



■0(e) 
(21) 



The first divergence, which is proportional to (q/B) 2 can be removed by wavefunction renormalization of the tree- 
level term ^B 2 in the free energy. This term is normally omitted since it is independent of the other parameters of 
the theory. The second divergent term, which is proportional to M 4 is the identical to the divergence that appears 
for B = 0. We can then add and subtract the term 



d d p 
J2nj ti 



vV + M 2 



M 4 



2(4tt) 2 



•0(e) 



(22) 



to Eq. (f2Tj) and take the limit d — > 3 in the difference. The divergence is now isolated in the integral on the left- 
hand-side of Eq. ([22]) . We set d = 3 here as well and and regulate it by imposing a sharp cutoff A in the usual way. 
The quark thermodynamic potential then becomes 



quark 



(M - m ) 5 



4G 

N c \- 

16tt 2 V 
/ 



Nc_ 
4tt 2 



AVA 2 + M 2 (2A 2 + M 2 ) + M 4 log 



A + VA 2 + M 2 



M 



(q f B) 2 [8C' (-!,£/) - 4(x 2 - x f )\ogx f + 2x 2 ] 



4 



E 

f,m,s 



\9fB\T r dp. 



2- 



2tt 



log 



l + 2$e 



-2/3E„ 



(23) 



The complete thermodynamic potential in a constant magnetic background is the sum of Eqs. (|13[) and (|23p and 
denoted by f2. The values of M and the Polyakov loop <E> are found by minimizing fl with respect to M and $, i. e. 
by solving the gap equations 



m = an 



(24) 



Using Eqs. (|23|) . we obtain 



(M - m ) _ Nc_ 
2ir 2 



2G 



a . iTH | A (2A + Af ) 2 
\/A 2 + M 2 



A + VA 2 + M 2 
M 



AT 



A 2 + M 2 + AVA 2 + Af 2 



167T' 



8C'(0,x/) - 2 (2a;/ - 1) \ogx f + 4x f 



+ 



Nc 
2tt 2 



E 



1_$2 



f,m,s 

24e~' 3Q 



-2/3 _E„ 



^ 2tt 

/,m,s 



E m 1 + 2$e-' 3B - + e- 2 ^- 



2tt 1 + 2<5>e-P E - 



e 



-20E„ 



(25) 



(26) 



We notice in particular that $ = is the only solution 
to Eq. flUD at T = and the PNJL model then reduces 
to the NJL model. 



III. NUMERICAL RESULTS 

The PNJL model has five different parameters, 
namely G, mo, and A in r2 qua rk and a and b in ri gaU gc- 
At T = 0, figauge = and so the PNJL model reduces to 
the NJL model. We can therefore determine the param- 
eters G, mo, and A separately. For A^ c = 3, one normally 
choses an ultraviolet cutoff A and tunes the parameters 
mo and G such that one reproduces the pion mass 
and the pion decay constant f„ in the vacuum. For 
N c = 2, we have no experiments toguide us and several 
different choices have been made [48|, HH, [6(| ■ We follow 
Ref. [48| that uses N c scaling arguments. Since the pion 
decay constant f„ is proportional to \/N c and the pion 
mass is proportional to N c , we simply rescale the three- 
color values by y| and |, respectively. This scaling 

gives U = 75.4 MeV and m^ = 93.3 MeV. Note that 
we in the following refer to the case where m^ =93.3 
MeV as the physical point. With an ultraviolet cutoff 
of A = 657 MeV, this gives G = 7.23 (GeV)~ 2 and 
m = 5.4 MeV at the physical point and G = 7.00 
(GeV)~ 2 and m = MeV in the chiral limit. The 



parameter a in f2 ga uge is related to the critical tempera- 
ture for deconfinement transition in pure-glue QCD and 
reads a — T c log24. In the pure gauge theory, T c is 
independent of the number of colors N c [6l| and so we 
can use the critical temperature for pure-glue QCD from 
lattice calculations with N c = 3, T c = 270 MeV. This 
yields a — 858.1 MeV. The parameter b can be tuned so 
that the chiral transition takes place at approximately 
the same temperature as the deconfinement transition. 
Note, however, that the Polyakov loop is strictly not an 
order parameter in the presence of dynamical fermions. 
It is like a crossover and the transition region is defined 
as a band in which $ varies. We define the transition 
region to be the temperatures where 0.4 < $ < 0.6 and 
the width of the band in the B—T plane tells of fast the 
crossover is. Nevertheless, we define a deconfinement 
temperature by the condition $ = | . This gives a curve 
in the B-T plane and acts a useful guide to the eye. 
The requirement $ = \ yields b = (210. 5) 3 (MeV) 3 . In 
Fig. [U we show the chiral condensate and the Polyakov 
loop in the chiral limit as a function of T/m v , where 
m,r = 140 MeV is the physical pion mass in the vac- 
uum. The chiral condesate (solid) line vanishes at the 
temperature at which $ = Thus for \qB\ = 0, the 
two transitions take place at the same temperature as 
explained above. For comparison, we also plot the chi- 
ral condensate in the NJL model as well as the Polyakov 
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loop in the pure-glue case, i.e. as derived from the po- 
tential rigaugc in Eq. (|13[) . The chiral transition in the 
chiral limit is second order for \qB\ = 0| in the NJL as 
well as in the PNJL model. 




TMr 

FIG. 1: Normalized chiral condensate and Polyakov loop as 
a function of T/m-n, where m-n = 140 MeV. 
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FIG. 2: The chiral condensate {'tptp) b normalized to the chi- 
ral condensate in the vacuum {ipip)o as a function of magnetic 
field \qB\ in the chiral limit (solid line) and at the physical 
point (dashed line) at zero temperature, m-* = 140 MeV. 



In Fig. [21 we show the chiral condensate (^ijj) b nor- 
malized to the chiral condensate in the vacuum as 
a function of \qB\ in the chiral limit and at the physical 
point in the vacuum i.e. T = fis = in the PNJL model. 
Note that here and in the following \qB\ is measured in 
units of the pion mass for N c = 3, i. e. = 140 MeV. 
Since the effects of the Polyakov vanish in the vacuum 
this is also the prediction of the NJL model. (ipip)B is 
a monotonically increasing function of \qB\ and the sys- 
tem exhibits magnetic catalysis. We notice that mag- 
netic catalysis is stronger in the chiral limit. 

We next consider the magnetic moment for a fermion 
of flavor /. In terms of the spin operator E Ml/ — 
Tih^l" - 7V), it is defined by (^/E^V/)- In the 
case of constant magnetic field in the z-direction, only 
E/ = E 12 is nonzero [63j. Using the properties of the 
7-matrices, it can be shown that only the lowest Landau 
level (LLL) contributes to the expectation value of E/ 
and reads 

$/E/^/>b = (^/>^ L \ (27) 

where the superscript indicates that we include only the 
lowest Landau level. We then define the polarization by 

fi r — — -p= 

{lpf1pf)B 

~ 1 77 / \ — ' ^> 

WfVf)B 

where the superscript HLL indicates that we have in- 
clude only the higher Landau levels. In Fig. [31 we show 



the polarization // = ^(/i u + fid) a t T = and in the 
chiral limit as a function of |<?.B|. As expected, the po- 
larization saturates for large magnetic fields to fi°° = 1. 
In this limit, the higher Landau levels effectively become 
very heavy (cf. Eq. ([T3])). they decouple and the LLL 
dominates the physics. In this limit all the fermions are 
in the LLL and their spin is pointing in the same direc- 
tion. 
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FIG. 3: The polarization ft as a function of \qB\ in the chiral 
limit and at zero temperature. m w = 140 MeV. 



In Fig. [4l we show the critical temperature for the 
chiral transition (solid line) and the critical tempera- 
ture for the deconfinement transition (dashed line) as 
functions of the magnetic field in the chiral limit for the 



G 



PNJL model. The band is denned by 0.4 < $ < 0.6 
and shows the transition region for the deconfinement 
transition. We also show the critical temperature in the 
NJL model for comparison. The parameter b in Eq. (fT5|) 
has been tuned such that the two transitions coincide for 
B = 0. We note that there is a splitting between the two 
transitions and that T c for the deconfinement is always 
lower than T c for the chiral transition and that the gap 
increases as a function of B. Thus, there should be a 
phase in which matter is deconfined and chiral symme- 
try is broken. The splitting was also observed in Ref. [29| 
where the authors coupled the Polyakov loop to linear 
sigma model with quarks with N c = 3 colors. More- 
over, while T c for the chiral transition increases by more 
than 20% from B = to \qB\ = 20ml, T d for the de- 
confinement transition is hardly affected. The width of 
the band is approximately 30 MeV. In contrast, the lat- 
tice simulations for N c — 2 reported in Ref. [55[ indicate 
that the critical temperature for deconfinement coincide 
with that of the chiral transition. Finally, we note that 
the determination of a and b has changes the critical 
temperature for the chiral transition dramatically. The 
increase of T c at B = is approximately 35 MeV and is 
fairly constant up to \qB\ = 20 MeV. In order to com- 
pare the chiral transition at finite magnetic field in the 
NJL and PNJL model, i. e. the effects of the Polyakov 
loop, there might be other ways of determining a and b. 
For example, one could force the deconfinement transi- 
tion and the chiral transition in the PNJL model to take 
place at the same temperature for B = and force it to 
coincide with the chiral transition in the NJL model as 
well. This way of determining the parameters in f2 ga ugc 
would not require the input from lattice simulations. 

In Fig.[5j we show the thermodynamic potential f2— fio 
divided by in the chiral limit for four different tem- 
peratures and \qB\ = 20mJ. For each temperature, we 
also give the value of the Polyakov loop. The critical 
temperature for the chiral transition is T c — 220 MeV 
and from the long-dashed line we see that transition is 
second order. Since the value of the Polyakov loop for 
T = 220 MeV is $ = 0.68, we conclude that the decon- 
finement transition has already taken place. 

In Fig.|6l we show the thermodynamic potential O — Oo 
divided by f* as a function of $ in the chiral limit for 
four different temperatures and \qB\ = 20m^.. For each 
temperature, we also give the value of the chiral conden- 
sate M. At T = 197 MeV, we find that the minimum of 
the effective potential occurs for $ = ^ which defines the 
critical temperature for the deconfinement transition in 
the chiral limit. At this temperature, M = 286.2 MeV 
and so we are still in the chirally broken phase. For 
T = 250 MeV, the chiral condensate is vanishing and so 
we are in the chirally symmetric deconfined phase. In 
the chiral limit, the chiral transition is always second 




FIG. 4: Critical temperature for the chiral transition (solid 
line) and critical temperature for the deconfinement transi- 
tion (dashed line) as functions of magnetic field \qB\ in the 
chiral limit for the PNJL model. The band is defined by 
0.4 < $ < 0.6. We also show the critical temperature for 
the chiral transition (thin line) in the NJL model, ro, = 140 
MeV. 
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FIG. 5: Thermodynamic potential — flo divided by as 
a function of M in the chiral limit for four different temper- 
atures and \qB\ = 2§m%, where m v = 140 MeV. See main 
text for details. 



order. 

At the physical point, the chiral transition is a 
crossover and there is no well-defined critical temper- 
ature T c at which the chiral condensate vanishes. How- 
ever, one can define a pseudo-critical temperature by the 
inflection point of the chiral condensate as a function of 
temperature. One can also define the transition region 
by the temperature range where M varies between 0.4 
and 0.6. This range depends on the magnetic field and 
gives rise to a band in the B-T plane. This is shown as 
the dark band in Fig. [7J In the same manner, we define 
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FIG. 6: Thermodynamic potential Q, — f2o divided by as 
a function of $ in the chiral limit for four different temper- 
atures and \qB\ = 20m^, where m n = 140 MeV. See main 
text for details. 
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FIG. 7: Transition region at the physical point as a function 
of \qB\. tUtt = 140 MeV. See main text for details. 



the crossover transition for the deconfinement transition 
by the the temperature range where 0.4 < $ < 0.6. This 
is shown as the light band in Fig. [7] Moreover, as a guide 
to the eye, we also show the lines where M = ^ and 
<J> = i, respectively. The curves indicate that the two 
transitions coincide for magnetic fields up to \qB\ « 5m \ 
and after that they split. The deconfinement transition 
is always taking place first. Finally, we note that the 
band of the chiral transition is much narrower than that 
of the deconfinement transition (approximately T = 10 
MeV versus approximately T — 30 MeV) and is so sig- 
nificantly faster. In Fig. El we show the normalized ther- 



modynamic potential fi//* as a function of M and $ in 
the chiral limit with \qB\ = 20m^ for different tempera- 
tures. From left to right T = 100 MeV, T = 200 MeV, 
T = 219 MeV, and T = 250 MeV. For T = 100 MeV, 
we are in the confined and the chirally broken phase, 
while for T = 200 MeV, the value of $ is slightly above 
\ ($ = | for T = 197 MeV). The minimum is still 
for nonzero M and so chiral symmetry is broken. For 
T = 219 MeV, M/M is still nonozero (T c = 220 MeV) 
and so we are just below the chiral transition. 



IV. CONCLUDING REMARKS 



In the present paper, we have considered two-color 
two-flavor QCD in a constant magnetic background us- 
ing the PNJL model. In the chiral limit, the chiral 
transition is always second order with mean-field crit- 
ical exponents. This is in agreement with universality 
arguments. Our result for T c for the chiral transition as 
a function of B for jV c = 2 is qualitatively the same as 
that for N c = 3, namely that it increases. For the PNJL 
model with lack of confinement and no dynamical gauge 
fields, we do not expect a qualitative different behav- 
ior. Whether the different results for T c for N c = 2 and 
N c = 3 obtained on the lattice is due to the different 
gauge groups and the very different physics in the vac- 
uum (diquarks rather baryons as bound states) or taking 
the continuum limit is an open question. For N c = 3, 
it seems that taking the continuum limit is essential to 



obtain the decrease of T c as function of B and so the 
only way to settle the issue for N c — 2 is to perform 
a careful study of the continuum limit in this case as 
well. Either way, given the lattice results for N c = 3 
it is clear that all the model calculations to date fail at 
temperatures T around T c and in particular does not 
incorporate that magnetic catalysis is turned into in- 
verse magnetic catalysis around the transition. This is 
independent of whether it is a mean-field calculation or 
goes beyond using e.g. functional renormalization group 
techniques [3^ 



18] 



1311 ] - However, in the recent paper 
the authors argue that fluctuations from neutral meson 
can exhibit magnetic inhibition, i. e. inverse magnetic 
catalysis if the magnetic field B is sufficiently strong. 
This may provide a mechanism that explains the recent 
lattice results. 



It is also interesting to note that the magnetic field 
hardly affects the critical temperature for deconfinc- 
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FIG. 8: The normalized thermodynamic potential fi//* as a function of M and $ for four different temperatures and 
\qB\ = 20m^. From left to right T = 100 MeV, T = 200, Mev, T = 219 MeV, and T = 250 MeV. See main text for details. 



merit. This is in line with the observation of .29j. Mor- 
ever, our results seem to indicate that the two transitions 
coincide at the physical point up to fairly large values of 
the magnetic field. \qB\ ~ 4m^.. 

Finally, we would like to comment on the role of quan- 
tum fluctuations and related renormalization issues. In a 
one-loop calculation of the effective potential it is possi- 
ble to separate the vacuum contributions from the ther- 
mal contributions. In some case, it therefore makes sense 
to investigate the role of the vacuum fluctuations. For 
example, in the QM, it is customary to treat the bosons 
at tree level and the fermions at the one-loop level. In 
this case it was shown in Ref. [64| for B = that the 
order of the phase transition depends whether the zero- 
temperature fluctuations are included or not; if they are, 



the chiral transition is second order and if they are not, 
it is first order. The same effect of the vacuum fluc- 
tuations were found in the entire \ib~T plane in strong 
magnetic fields in [32[ . In contrast to the QM model with 
quarks, this question does not make sense in the (P)NJL 
model. The reason is that chiral symmetry breaking in 
the (P)NJL is always a loop effect in contrast to the 
QM model where it is built into the tree-level poten- 
tial. In a similar manner, Ref. [65| finds that a crossover 
transition (for N c — 3) at the physical point remains 
a crossover at finite magnetic field B in an NJL model 
calculation. This is in contrast to Ref. [29j where it 
is found that strong magnetic fields turn the crossover 
into a first-order transition. In their work, the authors 
use the QM model and renormalize by subtracting the 
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fermionic vacuum fluctuations at B = 0. 
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